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CN , Abstract 

+^ , We measured the decay time of the scintination pulses produced by electron 

' O ■ and nuclear recoils in CaF2(Eu) by a new fitting method. In the recoil energy 

O 



region 5-30keVee, we found differences of the decay time between electron and 
nuclear recoil events. In the recoil energy region above 20 keVee, we found that 
the decay time is independent of the recoil energy. 
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1. Introduction 

> 

f~^ ' There is convincing evidence that most of the matter in our Galaxy must be 



in 



dark matter. One of the most popular candidates for dark matter is Weakly In- 
teracting Massive Particles (WIMPs). WIMPs are thought to be non-baryonic 



f— ^ ■ particles, and the most plausible candidates for them are the Lightest Super- 

^"^ ' symmetric Particles (LSPs) and the Lightest Kaluza-Klein Particles (LKPs)[l|. 

WIMPs can be searched for directly, as they interact with atomic nuclei in 
K> , detectors. Direct detection relies on one of two modes of interaction with target 

jrt ■ nuclei. The first mode is called spin-independent (SI) coupling. It describes 

coherent interaction with the entire nuclear mass. Therefore the SI part of 
the WIMP-nucleus cross section is large, if the target nuclei have the large 
mass number. For example, xenon and iodine are the favorable nuclei to detect 
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WIMPs for the SI interaction. The second mode is cahed spin-dependent (SD) 
couphng. It describes interaction of WIMP with the spin-content of the nucleus. 
Hence the SD part of the WIMP-nucleus cross section is large, if the spin-content 
of the target nuclei is large. For example, ^^F is one of the most favorable nuclei 
to detect WIMPs for the SD interaction because of its large nuclear spin. 

Several direct WIMP searches using^^F-based detectors, such as bolometers M^. 
3|, a bubble chamber [jl, scintillators [5|, la, |7|,|8|,|9|, a superheated droplet detector 
and so on, have already been performed. In our group, a WIMP search exper- 



io| 



iment using a CaF2(Eu) scintillator is carried out at Kamioka observatory in 
2005[5|- CaF2(Eu) is one of the most suitable material among ^^F-based scin- 
tillator because of its high light-output. 

In other direct WIMP searches, pulse shape discrimination (PSD) technique 
is used in order to statistically discriminate nuclear recoil signal events from 
electron recoil background events. A nuclear recoil event is different from an 
electron recoil event in decay time of their pulse shapes. In the PSD technique, 
a distribution of the decay time constant is utilized for the statistical discrimi- 
nation. The PSD feature in CaF2(Eu) was reported in the MeV region[9l lll| . 
The pulse shape of CaF2 (Eu) scintillator in lower energy region relevant to the 
dark matter WIMP search was measured previously by Tovey et al. [12| . They 
reported that CaF2(Eu) has no PSD capability in the 10-30 keVee region. How- 
ever, we recently measured the pulse shape of CaF2(Eu) and analysed our data 
by a more careful statistical method. In this paper, we report on the new mea- 
surement of CaF2(Eu) pulse shape and the analysis. We found a difference in 
the pulse shape of CaF2(Eu) between nuclear and electron recoil events. 

2. Experimental methods 

2.1. Measurement system 

The scintillator sample consisted of an unencapsulated, double-ended, 20- 
mm-diameter, 20-mm-long crystal of CaF2(Eu) with europium doping content 
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Figure 1: Schematic view of tiie neutron detectors. The neutron events is selected by the 
recoil energy in CaF2(Eu) scintillator and the TOP between the plastic scintillator and the 
liquid scintillator. 



of 0.5%. The crystal was glued on two 2-incli photomultiplier tubes (Hamamatsu 
R329-02). 

^^^Cf was used as a neutron source together with three detectors as shown 
in Fig.[T] Nuclear recoil candidates were selected by looking for coincidences 
between events in three detectors. The first one was the target CaF2(Eu). The 
second was a plastic scintillator, which detected prompt 7 rays from fission 
decays of ^^^Cf. The last was a liquid scintillator (Saint-Gobain, BC-501A), 
which caught scattered neutrons from the target CaF2(Eu). 

For the latter two detectors, Hamamatsu R329-02 photomultiplier tubes 
were also used. Signals from all the photomultiplier tubes were sent to fast 
discriminators with sufficiently low thresholds. A special care was taken to 
the two output signals of the CaF2(Eu) scintillator so as to suppress multiple 
trailing pulses due to its long decay time constant. A coincidence of these two 
discriminator outputs were taken to define the CaF2(Eu) scintillator hit with a 
coincidence width of 160 ns. Coincidence widths of the plastic scintillator and 
CaF2(Eu)scintillator, and the plastic scintillator and the liquid scintillator were 
set to 195 ns and 400 ns, respectively. 



The photomultiplier outputs of the CaF2(Eu) and the hquid scintniators 
and the discriminator output of the plastic scintihator were recorded with a 
Tectronix TDS3034B digital oscilloscope with sampling frequency of 2.5 GHz 
and 300 MHz bandwidth by a trigger of the coincidence of three detectors. The 
pulse shape of the CaF2(Eu) scintillator and the timing of all the detectors are 
thus recorded and sent to a computer. 

The estimated solid angles of the CaF2(Eu) scintillator as seen from the 
^^^Cf source and of the liquid scintillator from the CaF2(Eu) scintillator are 
1.2 X 10^"^ sr and 1.1 x 10^^ sr, respectively. The scattering angle spanned 
by the three detectors was about 90 degrees. Between the plastic and liquid 
scintillators, polyethylene (PE) blocks (35 cm long) and Pb blocks (20 cm long) 
are placed in order to shield the liquid scintillator from neutrons coming directly 
from 252 Qf 

2.2. Selection of nuclear recoils 

We distinguished nuclear recoil events from electron recoil events using the 
time-of-flight (TOE) of the particle. The time-of-flight of the neutron from the 
252 Cf source to the liquid scintillator is 

tsL — tsT + tTL = — , H , (1) 

Here tsT and ixL are the times-of- flight of the neutron from the 252 Cf source to 
the target scintillator and from the target scintillator to the liquid scintillator, 
and LsT and Ltl are the respective lengths. In this measurement, Lst and 
LxL are 70 cm and 15 cm. Ei and E2 are the neutron energies before and after 
the collision with the nucleus in the target, and m^ is the neutron mass. Then, 
the energy of the recoil nucleus is 

2?Ti,i 



Er = Ei-E2^EiX ,,^ ;•'" X [M + mn sin2 - cos 9 J M"^ - ml sm^ ( 

^n) V 



(2) 

when the scattering is elastic. Here M is the recoil nucleus mass, and 6 is the 
neutron scattering angle. In this measurement, since the recoil angle is around 
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Figure 2: Scintillation energy and TOF distribution by density plot. The curves are the 
theoretical values of the scattering events for the ^^F and ''"Ca with quenching factors of 0.05 
and 0.15. The inside of the rectangle is selected as neutron events. 



90 degrees, the recoil energy is simply 



E, = E, 



E. 



E^ X 



2m„ 



(3) 



n 



M + rrin 

Inelastic cross section is negligibly small in the neutron energy range [l3| relevant 
to the present measurement and therefore ignored in the analysis. 

Density plots of the electron equivalent energy deposited in CaF2(Eu) versus 
TOF tsL is shown in Fig.[2l Theoretical expectation curves are also show in the 
Fig. [5] for the ^^F and "^^Ca with quenching factors of 0.05 and 0.15. On the 
other hand, Compton scattering events in the CaF2(Eu) scintillator caused by 
gamma rays emitted by the ^^^Cf source are expected to gather in the region 
where TOF is ns. Fig.[2]shows clear separation of the nuclear recoil events and 
the Compton scattering events. 

To pickup the ^^F and "^^Ca recoil events, we put a rectangle where nuclear 
recoil events are expected. Consequently, we obtained enough pulse shapes for 
the PSD analysis. 



2.3. New method of PSD 

A pulse shape of scintillation is formed by a train of small single-photoelectron 
pulses. The arrival time relative to the pulse start time of each pulse follows 
an exponential distribution. Hence, the scintillation pulse shape has a char- 
acteristic time constant, which is different between electron and nuclear recoil 
events. 

In many researches, integrated pulse is used for analysis, and an exponential 
rise function Ax (1 — exp[— (t — to)/T]) is fitted to the data by the least-squares 
method|14|. This analysis is statistically not strictly correct because integrated 
data points are not independent of each other. Therefore, the data themselves 
before integration must be used for the fitting analysis. 

However, because the scintillation of CaF2(Eu) is not bright enough and 
has large decay constant, it is difficult to fit an exponential fall function A x 
exp[— (t — to)/T] to all the data points by the least-squares method. So we 
divided a train of data points of the oscilloscope into 40 ns bins, and fitted the 
estimate number of photoelectrons in each bin calculated by an exponential fall 
function to the data by the maximum-likelihood method via two parameters, 
the estimate value of the first bin and the decay time constant. 

The number of photoelectrons in each bin follows a Poisson distribution. So 
the probability distribution function f{xi, fit) of the observed charge Xi in i-th 
bin is 

f{xi, ^JLi) = ^ [Poi(n, ^li) ■ Gau(a;i, Qo + nQi.al + nal)] , (4) 

n 

where /^i is the estimate value of the number of photoelectrons in the i-th bin. 
Poi(n, /i) is a Poisson distribution function of n with the mean value of /i, and 
Gau(a;, Q, cr^) is a Gaussian distribution of x with the mean value of Q and 
dispersion a^ . Qq, (Tqi Qi and cti are the parameters for the conversion of 
the number of photoelectrons to the charge recorded by the oscilloscope. Qo 
and (To a-re the pedestal and its width, and Qi and cti are one photoelectron 
response and its standard deviation. These four parameters are characteristic 
of PMTs and they were measured by a single photoelectron spectra analysis 



independently [15|. Tab. [T] shows the typical values of Qq, o-q, Qi and cti as 
measured in electron equivalent deposit energy. The Qi, ctq and cti values 
are stable within ±4.5%, ±16% and ±23% in measurement after over a year, 
respectively. 

Table 1: Values of Qg, ctq, Qi and ai as measured in electron equivalent deposit energy. 

Parameter PMT ch.l [keV] PMT ch.2 [keV] 

go (1.1 ± 0.3) X 10-3 (2.2 ± 0.3) X 10-3 

ao (1.56 ± 0.04) X 10-2 (1.83 ± 0.04) x 10"^ 

Qi (3.15 ± 0.07) X 10-1 (3.06 ± 0.07) x lO^^ 

(1.70 ± 0.08) X 10-1 (1.69 ± 0.08) x 10"! 



(Tl 



From this, the likelihood function L{fiQ,T) is 

PMTch.1,2 i ^ \ ''' 

where t is the bin width. The two parameters, /iq (the estimate value of the 
photoelectron number in the first bin) and t (the decay time constant), are 
evaluated when they maximize the likelihood function L(fio,T). 

3. Result and discussion 

A sample of the result of the maximum likelihood fitting to the experimental 
data is shown in Fig. [3] for 7 ray events with ^''^Am and ^^^Cs, and for selected 
neutron events with ^^^Cf. A difference is seen in the decay time distribution 
between electron and nuclear recoil events. 

In order to see the difference more clearly, we divided the events into elec- 
tron equivalent energy bins of 5 keV width, and evaluated the mean value of 
the decay constant in each bin. Exceptionally, we put regions of photoelectric 
absorption peaks, 13.9keVee (^^lAm), 59.5keVee (^^lAm) and 32.2-36.4keVee 
(13'i'Cs), into wider single bins. The distribution of the decay time constant r 
can be approximated by a Gaussian in ln(T)|l6l|: 
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Figure 3: The scatter plot of the experimental data. The decay constant is evaluated by the 
maximum likelihood method. ('Sec. |2.3T l 




gamma from Am241 ^ 
gamma from Cs137 
neutron from Cf252 ^ 
UK: gamma > 
UK: neutron 



-^trr 



30 40 50 

Energy [keVee] 



Figure 4: The decay time constant versus the energy. The present result is shown by open 
circle, filled circle and open triangle signs. The results in Ref. |12l | are open and filled squares. 



We evaluate three parameters, r^ (the mean value of the exponent), N^ (the 
number of events) and w (a width parameter), in each bin by the least-squares 
method. Fig. 2] shows the evaluated mean decay time constants t^ as a function 
of energy. There are marked differences of the decay time between electron and 
nuclear recoil events. 

The results in Ref. \l2\ are also shown in Fig.|31 In their results, the differ- 
ences are not seen. Moreover, the decay constants rise rapidly with energy from 
~ 660 ns to ^ 850 ns. In the present work, a similar dependence exists in the 
low energy region, but it is weak. 

Strictly speaking, because the pulse shape is not a single exponential curve, 
but the sum of a few exponential curves, the value of the decay time constant is 
dependent on the way of the analysis. So, the values of different measurements 
cannot be compared directly. However, these properties, the difference between 
the nuclear recoil and the electron recoil and the relation between the decay 
constant and the recoil energy, are independent of the way of the analysis. 

4. Conclusion 

We developed a new method to evaluate the decay time constant of scintilla- 
tion and applied it to the CaF2(Eu) scintillator. We found the difference in the 
scintillation decay constants between electron and nuclear recoil events (Fig.S]). 
We also concluded that the decay time constant of CaF2 (Eu) is constant in the 
recoil energy region above 20keVee. 
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